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Abstract 

We consider the conditions for the decay products of perturbative inflaton 
decay to thermalize. The importance of considering the full spectrum of inflaton 
decay products in the thermalization process is emphasized. It is shown that the 
delay between the end of inflaton decay and thermalization allows the thermal 
gravitino upper bound on the reheating temperature to be raised from 10 8 GeV 
to as much as 10 12 GeV in realistic inflation models. Requiring that thermaliza- 
tion occurs before nucleosynthesis imposes an upper bound on the inflaton mass 
as a function of the reheating temperature, m,s ^ 10 10 (Tr/1 GeV) 7 / 9 GeV. It 
is also shown that even in realistic inflation models with relatively large reheat- 
ing temperatures, it is non-trivial to have thermalization before the electroweak 
phase transition temperature. Therefore the thermal history of the Universe is 
very sensitive to details of the inflation model. 
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1 Introduction 



Reheating is a fundamental process in early Universe cosmology [[jj , in which the energy 
density in a coherently oscillating inflaton field is converted to thermalized relativistic 
particles. Originally it was believed that this occurred simply by the perturbative de- 
cay of the individual scalar particles in the corresponding Bose condensate of inflatons 
]]], 0, but in recent times it has become clear that the process can be considerably 
more complicated, with non-perturbative processes such as parametric resonance || 
and quantum creation of fermions [fj. [5|. |(| playing a significant role. Nevertheless, in 
many models, even if there is an early stage of preheating, the latter stage of reheat- 
ing is dominated by the perturbative decay of the remaining inflaton energy density. 
An important issue is then the thermalization of the inflaton decay products. Since 
the inflaton can be a very massive particle, as heavy as 10 15 GeV in some typical 
inflation models ||, it is not obvious that its highly energetic decay products will 
thermalize rapidly. In this paper we will consider the conditions under which complete 
thermalization of the inflaton decay products occurs^]. We will see that the results are 
particularly important in SUSY inflation models, allowing the thermal gravitino upper 



bound on the reheating temperature |L0], [LI] to be substantially increased. In addition, 
we will show that thermalization can occur at low temperatures even in realistic infla- 
tion models, as low as the electroweak phase transition temperature or less. We will 
also calculate the upper bounds on the inflaton mass following from the requirement 
of thermalization before nucleosynthesis |TT|| . 



2 Thermalization of Inflaton Decay Products 

There are two distinct processes involved in thermalization of the decay products, 
which we shall refer to as " self-thermalization" and "catalysed thermalization". 

In order to avoid confusion, we first define what we mean by the "thermaliza- 
tion temperature" and the "reheating temperature". We will define radiation as a 
1 An earlier discussion of the thermalization of inflaton decay products is given in fl. 



1 



background of relativistic particles. The thermalization temperature, T t h, refers to 
the temperature of the radiation when the relativistic particles can scatter rapidly 
enough relative to the expansion rate of the Universe to come into thermal equilib- 
rium. The reheating temperature, Tr, is defined as the temperature the radiation 
would have at the time when the Universe becomes radiation dominated if it were 
in thermal equilibrium. This is the conventional reheating temperature of inflation 
models (T R « (M P iT d ) l l 2 , where T d is the inflaton decay rate), which usually assume 
that thermalization of the inflaton decay products is instantaneous, 
(i) Self-Thermalization 

The energy density of the decaying inflaton field, S, is given by 



where a is the scale factor. (We are assuming here that a single decaying inflaton field is 
the source of the thermal energy.) Thus most of the energy density in the inflaton field 
decays when H « (where H = 2/3t is the expansion rate during inflaton matter 
domination), just before the Universe becomes dominated by relativistic particles. 
Therefore the apparent condition for the inflaton decay products to thermalize is that 
these decay products should thermalize by scattering from each other. (This is the 
condition considered in ||.) We refer to this process, the thermalization of decay 
products produced during an interval St ~ H^ 1 by scattering from each other, as self- 
thermalization. However, we will see later that this condition for the thermalization 
of the radiation background is incorrect; there are also much lower energy particles in 
the spectrum of decay products, coming from the red-shifted decay products of earlier 
inflaton decays, which play a crucial role in the thermalization process. 

We first derive an upper bound on the inflaton mass from self-thermalization of 
the decay products produced during St ~ H -1 at the end of inflaton decay. The 
initial energy of the decay products will be of the order of the inflaton mass, ms- The 
condition for the thermalization of these decay products by scattering from each other 
is then 




(1) 



An(H)a sc (H) 



N SC H , 



(2) 
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where N sc is the number of scatterings required to fully thermalize the energy; typically 
N sc < 10. An is the number of decay products at H which were produced in a time 
St ~ Hft 1 at H R , where H R is the expansion rate when the Universe becomes radiation 
dominated. This is given by 

, ■ >, 

a H 



An(H) » f^Y r d H^n s , (3) 
V an / 



where ris ~ ps(H R )/m s is the number of inflatons remaining in the condensate at H R , 
an is the scale factor at H and Yd is the inflaton decay rate, given by 

where T R is the conventional reheating temperature and g(T) is the number of rela- 
tivists degrees of freedom. The scattering rate at H for relativistic particles of initial 
energy E fa ms at H R is given by 

a 2 ( a \ 2 a 2 



- \aj m 2 ' (5) 

where a = g 2 j^n corresponds to the gauge or Yukawa couplings. (For now we will 
consider massless decay products.) This assumes that 2 — > 2 particle scattering pro- 
cesses can produce final state particles which subsequently rapidly decay, so increasing 
the number density and decreasing the average energy of the particles in the ensemble; 
otherwise we should consider processes such as 2 — > 4 particles, with a correspond- 
ingly smaller a. Thus the condition for complete self-thermalization at a temperature 
T <T R is 

< (T R \^ f?>M Pl k T a 2 \ 1,Z ^ 2/3 



Numerically we find 

m „„ = 2, x 10' aV (^) " 3 (j^) GeV , (7) 

where we have used kr ~ 17 (g(T) fa 100). 
(ii) Catalysed Thermalization 
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The naive approach to thermalization considers only the self-thermalization of the 
decay products produced at the end of reheating. In fact, there will also be red-shifted 
decay products from earlier inflaton decays. If these have red-shifted sufficiently, their 
scattering rate (oc E~ 2 ) can become large and so they can self-thermalize, transferring 
their energy density from a small number of high energy particles to a larger number 
of low energy particles. These can then act as targets for higher energy particles 
in the energy spectrum to scatter from and thermalize, with the process continuing 
until all the decay products are thermalized. We refer to this process as "catalysed 
thermalization" . The conditions for catalysed thermalization to occur are then that (i) 
there are particles in the energy spectrum of decay products of sufficiently low energy 
as to be able to self-thermalize and so provide a "seed" for catalysed thermalization 
and (ii) that catalysed thermalization can then thermalize the whole spectrum of decay 
products in a time < H~ l . 

In order to discuss catalysed thermalization, we need the spectrum of decay prod- 
ucts, dn(E,T)/dE, at T. Inflaton decay during a time 5U m H~ l at Hi contributes a 
number density at H given by 

dn(H,H t )~( a ^)\ d Hr P -^. (8) 
V an J m s 

The energy of the decay products red-shifts to 

"# R \ 2/3 /#\ 1/2 



E 



a Hi \ (H R y*(Hy* 



for H < H R . In the following, it will be sufficient to consider the spectrum at H < H R , 
since the weakest bounds generally correspond to both the largest red-shift of the decay 
products and the smallest value of H . Thus we find 

\T R J m s \m s J 

The change in energy at H of the decay products produced at Hi in a time 5ti ss H~ x 
is 5E w 2E/3. Thus 

dn ^3(T\W PS (H R ) EW 

dE~2\T R ) m s m 3/2> W 
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for H < Hr. An important point in what follows is that this spectrum has a low 
energy cut-off, at E min , corresponding the the inflaton decay products produced at 
the earliest time, immediately after the end of inflation at H = Hi, 

fH n \ 2 ' 3 /H \ 1/2 / T \ I k TE rt\ m , v 

The condition for a self-thermalized seed to exist at H < Hr is then that, for some 
energy E c > E min , self-thermalization of the decay products can occur for all E up to 

E c , 

dn a 2 ^ , . , 

— - > N SC H , WE < E c . (13) 

Using Eq. (pTTf) , we find that E c is given by 

9^(T R )« 4 M^T| 

320m|iV s 2 c T ' 1 J 

A self-thermalized seed will therefore exist if E m i n is less than E c at the smallest value 
of T, which imposes an upper bound on the inflaton mass 



< / t^ 1 ' 3 ( 9^(r H )ji4(v \ 1/6 4/9 1/9 _ 

mS ~W i 320^ J T « if '= m «-" (15) 
where T <Tr. Numerically we find 

- U * ^ W (l^f QeV . (16) 

If this is satisfied, then catalysed thermalization can thermalize the spectrum for all 
E such that 



_E cm = y/EEth is the centre of mass energy for scattering between an inflaton decay 
product of energy E and a thermalized particle of energy E t h ~ Pi/ 4 , where pe ~ 
Edn/dE is the energy density in particles of energy i? to 2E, and we have taken the 
number density of thermalized particles to be « p^ 4 . The condition for catalysed 
thermalization is then that 



dn . /N.,H 
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Since dn/dE increases with E, the condition for catalysed thermalization will be sat- 
isfied for all E > E c so long as it is satisfied at E c , which requires that E t h(E) < E c 
at E — E c . This is true so long as 

ms < 0.6a^N-^g(T R )^ f^f^j = m cat . (19) 

Numerically we find 

m cat « 1.7 x 10 13 a 6 / 5 iV- 3 / 5 (-g^) (i^)' 75 GeV . (20) 

As this is a much weaker upper bound than m see d, the upper bound from catalysed 
thermalization is generally given by mg ^ m S eed- The upper bound from catalysed 
thermalization is considerably weaker than the naive upper bound based on self- 
thermalization, Eq. (|7|), typically by two to three orders of magnitude. 

3 Consequences for Thermal Gravitinos, Nucleosyn- 
thesis and the Electroweak Transition 

So far we have not considered a specific inflation model, so our results are valid for both 
SUSY and non-SUSY models. In general, there is a lower bound on the thermalization 
temperature from nucleosynthesis, T t h ^ 1 MeV jTT]. In addition, in SUSY inflation 
models there is an upper bound from requiring that gravitinos are not produced ther- 
mally, T t h ~ 10 8-9 GeV for gravitino masses in the range 100 — 500 GeV |T^, O. 
Usually it is assumed that the inflaton decay products thermalize instantaneously, so 
that T t h is identified with Tr. However, this depends on the inflaton mass. Perhaps 
the most interesting consequence of this is that the thermal gravitino upper bound on 
T R can be considerably relaxed in realistic SUSY inflation models. To see this, we note 
that thermal gravitinos can only be generated at T ^ T t h- Thus the thermal gravitino 
upper bound should be T t h ^ 10 8-9 GeV. If ms ^ rn see( i when T m 10 8-9 GeV, then 
thermalization will occur safely below the thermal gravitino upper bound for the cor- 
responding value of T R . In Table 1 we give values of m seed as a function of T R for the 
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case T = 10 8 GeV. (The values of m see d for T = 10 9 GeV are given by multiplying 
the values in Table 1 by 0.46.) From this we see that for inflaton masses in the range 
10 15-16 GeV, as would be expected, for example, in D-term inflation models^, the up- 
per bound on the reheating temperature Tr is 10 10-12 GeV. Given the importance of 
the thermal gravitino upper bound as a constraint on inflation models, this weakening 
of the upper bound on Tr is significant. 



Table 1. Inflaton Mass Lower Bounds from Thermal Gravitino 

Non-production. 



Tr 


mW(« 2/3 iV- 1/3 


VlO 13 GeV J > 


10 s GeV 


3.2 x 10 13 


GeV 


10 9 GeV 


1.9 x 10 14 


GeV 


10 10 GeV 


1.1 x 10 15 


GeV 


10 11 GeV 


6.9 x 10 15 


GeV 


10 12 GeV 


4.1 x 10 16 


GeV 


10 13 GeV 


2.5 x 10 17 


GeV 


10 14 GeV 


1.5 x 10 18 


GeV 


10 15 GeV 


8.9 x 10 18 


GeV 



The nucleosynthesis lower bound on T t h imposes upper bounds on the inflaton 
mass. These bounds are relatively weak for large Tr, but for smaller values of Tr 
they can be significant. In the context of SUSY models there have been some recent 
motivations for considering low reheating temperatures. One is from Affleck-Dine 



baryogenesis p3| , |T4| , |15j . For the lowest dimension R-parity conserving flat directions 
of the MSSM scalar potential, those with dimension d = 4 and 6 (where the dimension 
refers to the non-renormalizable superpotential terms responsible for lifting the flat 
directions |L4|, |l|, [16[|), the observed baryon asymmetry requires that the reheating 
temperature is approximately 10 7 GeV and 1 GeV respectively |15|, [0J. So Tr « 



2 The inflaton mass in D-term inflation models is given by ms = A£, where A is the Yukawa 
coupling of the inflation sector fields and the microwave backgoriund implies that («7x 10 15 GeV 
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1 GeV is one favoured possibility if the baryon asymmetry originates via the Affleck- 
Dine mechanism. Another motivation for low reheating temperatures is the possibility 
that large non-thermal gravitino densities are created by the oscillating inflaton field 
at the end of inflation ||. Although the resulting upper bound on the reheating 
temperature is sensitive to the details of the inflation model, there are indications 
that, for mass scales typical of inflation models, the upper bound is likely to bef] 
Tr <10 3 GeV §. 

Table 2. Inflaton Mass Upper Bounds vs. Thermalization Temperature 



Tr 


Tth 


m self /a 2/3 


m see ,/(« 2/3 iV- 1/3 (To#kv) V9 ) 


1 GeV 
1 GeV 


1 MeV 
1 GeV 


2.9 x 10 7 GeV 
2.9 x 10 6 GeV 


9.2 x 10 10 GeV 
9.2 x 10 9 GeV 


10 3 GeV 
10 3 GeV 
10 3 GeV 


1 MeV 

10 2 GeV 

10 3 GeV 


2.9 x 10 10 GeV 
6.3 x 10 8 GeV 
2.9 x 10 8 GeV 


2.0 x 10 13 GeV 
4.3 x 10 11 GeV 
2.0 x 10 11 GeV 


10 8 GeV 
10 8 GeV 
10 8 GeV 


1 MeV 
10 2 GeV 
10 8 GeV 


2.9 x 10 15 GeV 
6.3 x 10 13 GeV 
6.3 x 10 11 GeV 


1.5 x 10 17 GeV 
3.2 x 10 15 GeV 
3.2 x 10 13 GeV 



In Table 2 we give the upper bound on the inflaton mass as a function of the 
thermalization temperature for Tr = 1 GeV, 10 3 GeV, and 10 8 GeV. For Tr « 
1 GeV the nucleosynthesis upper bound from catalysed thermalization implies that 
ms ^ 10 10 ~ n GeV. Given that the inflaton mass scale in SUSY inflation models can 
naturally be ms ~ 10 15 ~ 16 GeV, this can impose a significant constraint on inflation 
models compatible with d = 6 AD baryogenesis. We also give the upper bound for 
the case T = Tr, corresponding to the case where the inflaton decay products instan- 
taneously thermalize. This requires that ms ^ 10 9 ~ 10 GeV for Tr « 1 GeV. These 
bounds are much weaker than would be expected from naive self-thermalization; com- 
paring m see d with m se if shows that self-thermalization would impose an upper bound 
3 This estimate is based on a model with a single chiral superfield, so that the spin-1/2 components 
of the gravitino are effectively the inflatino || . It remains to be seen whether this remains true in 
more realistic models. 
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smaller by a factor of more than 10 3 . For the case T R fa 10 3 GeV, the nucleosynthe- 
sis upper bound requires that mj < 10 12 ~ 13 GeV, which is less than would typically 
be expected in many inflation models, but which could nevertheless be satisfied with 
some moderately small couplings. For the case T R fa 10 3 GeV we have also calculated 
bounds for T fa 10 2 GeV, corresponding to thermalization before the electroweak 
phase transition temperature, which is necessary for the existence of an electroweak 
phase transition. We see that this imposes quite a strong upper bound on the infla- 
ton mass, ms ^ 10 11 GeV. Finally, instantaneous thermalization for T R ~ 10 3 GeV 
requires that m s < 10 10 " 11 GeV. For the case where T R is of the order of the con- 
ventional thermal gravitino upper bound, T R fa 10 8 GeV, the nucleosynthesis bound 
requires that ms 1z 10 16 ~ 17 GeV, which is easily (although not necessarily trivially) sat- 
isfied in inflation models such as D-term inflation. The electroweak transition bound 
requires that m s < 10 14 ~ 15 GeV. Again, although this can be satisfied in many in- 
flation models, it is nevertheless of the same order of magnitude as the inflaton mass 
expected in D-term inflation models, for example. Therefore it is non-trivial to have 
thermalization before the electroweak phase transition temperature, even in realistic 
inflation models with relatively high reheating temperatures! Finally, instantaneous 
thermalization for T R fa 10 8 GeV occurs only if ms "Z 10 12-13 GeV. 

This shows that the assumption of a thermalized background following inflation, 
and in particular the existence of an electroweak phase transition, strongly depends 
upon the inflaton mass. In many inflation models it is likely that the relativistic parti- 
cles of the radiation background will only thermalize at a relatively low temperature. 
Although for larger reheating temperatures it is quite easy to have thermalization be- 
fore nucleosynthesis, it is not so clear that thermalization will occur before the epoch 
of the electroweak phase transition. For lower reheating temperatures, thermalization 
before nucleosynthesis can impose significant bounds on the inflaton mass and so on 
the model of inflation. 

Finally, we comment on some assumptions made in this discussion. We have as- 
sumed throughout that the inflation decay products are massless. Of course, only the 
photons and possibly the gluons and neutrinos can be treated as massless throughout. 
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However, this is sufficient for our discussion of thermalization, since once these parti- 
cles thermalize they will serve as a thermalized background which can thermalize all 
the other particles. In addition, we have also assumed that the only source of radi- 
ation is the single decaying inflaton field. In fact, shortly after the end of inflation, 
there may be other sources of radiation which could contribute to the thermalized 
seed leading to catalysed thermalization. This could provide a much larger number of 
seed particles, relaxing the upper bound from m see d- For example, in D-term inflation 
models, in addition to the inflaton there are the fields _ responsible for hybrid 
inflation 0. The rapid decay of the ip- field (which has a mass typically of the same 
magintude as the inflaton mass |12] ) releases roughly the same energy in a time 
5t ~ Hj l as that stored in the inflaton field at the end of inflation. Therefore the 
number of lowest energy particles in the spectrum, corresponding to those produced 
at if ~ Hi, is enhanced by a factor Hi/ Hr. The result is that the m see d upper bound 
is increased by a factor (Hj/ Hr) 1 / 3 , which is greater than 10 5 for Hi « 10 13 GeV and 
Tr ~Z 10 8 GeV. Thus the bounds on the inflaton mass in hybrid inflation models are 
likely to be much weaker than in models with a single inflaton field. 

4 Conclusions 

We have considered the constraints following from the requirement that the relativistic 
decay products from inflaton decay thermalize. We have shown that the low energy 
decay products from inflaton decays occuring shortly after the end of inflation play a 
vital role in the thermalization of the whole energy spectrum, a process we refer to 
as catalysed thermalization. For the case where a single decaying inflaton field is the 
source of the thermal energy, requiring that the decay products do not thermalize be- 
fore the temperature of the thermal gravitino upper bound allows the upper bound on 
the "reheating temperature" in SUSY inflation models to be increased from 10 8 GeV 
to 10 10-12 GeV in realistic inflation models. Requiring that thermalization occurs be- 
fore to onset of nucleosynthesis can impose tight upper bounds on the inflaton mass 
for low reheating temperatures, such as may be required by Affleck-Dine baryogenesis 
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or non-thermal gravitino production. In addition, even in realistic inflation models 
with a relatively high reheating temperature, it is quite possible that the relativistic 
background will not thermalize before the temperature of the electroweak phase tran- 
sition. In hybrid inflation models, such as D-term inflation, the thermal energy due 
to the decay of the hybrid inflation fields at the end of inflation can provide a seed of 
thermalized particles which makes subsequent thermalization of the inflaton's decay 
products much more efficient, relaxing the upper bounds on the inflaton mass. The 
lesson from all this is that in many inflation models the thermal history of the Universe 
is likely to be quite different from that naively assumed on the basis of instantaneous 
thermalization of the inflaton's decay products. 
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